Little-Parks Oscillations in Hybrid Ferromagnet-Superconductor Systems 
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On the basis of of linearized Usadel equations we consider superconductivity nucleation in multi- 
ply connected mesoscopic superconductor/ferromagnet hybrids such as thin- walled superconducting 
cylinders placed in electrical contact with a ferromagnetic metal. We study the interplay between 
the oscillations of T c due to the Little-Parks effect and the oscillations due to the exchange field. 
We demonstrate that the exchange field provokes the switching between the superconducting states 
with different vorticities and this may result in the increase the critical temperature of the supercon- 
ducting transition in the magnetic field. Moreover we analyse the influence of the S/F transparency 
on the realisation of the states with higher vorticities. 

PACS numbers: 74.25.Dw, 74.45. +c, 74.78.Na 



I. INTRODUCTION 

Little-Parks effect^ i.e., the oscillations of the criti- 
cal temperature T c of multiply-connected superconduct- 
ing samples in an applied magnetic field H, is one of 
the striking phenomena demonstrating coherent nature 
of the superconducting state. Nevertheless, such oscilla- 
tory phase transition line T C {H) was shown to be inherent 
to simply connected mesoscopic samples with the lateral 
size of the order of several coherence lengths {&2£ and to 
hybrid ferromagnet (F)-superconductor (S) systems with 
magnetic dots or domains^i^, that create a "magnetic 
template" for nucleation of the superconducting order pa- 
rameter. These T C (H) oscillations reflect the switching 
between vortex states characterized by different winding 
numbers and are associated with the orbital effect. 

Another mechanism of switching between the super- 
conducting states with different vorticities in multiply 
connected hybrid S/F structures such as thin- walled su- 
perconducting shell placed in electrical contact with a 
ferromagnetic cylinder was suggested recently ini£. This 
mechanism is caused by the exchange interaction and 
associated with the damped-oscillatory behavior of the 
Cooper pair wave function in a ferromagneliii. It was 
shown that under certain conditions the exchange in- 
teraction can stimulate the superconducting states with 
a nonzero vorticity in the absence of external magnetic 
field. The interplay between the exchange and orbital 
effects may result in a subsequent switching between the 
states with different vorticities, as the F core radius in- 
creases. An obvious consequence of these transitions be- 
tween the states with different L should be a nonmono- 
tonic dependence of the critical temperature T c on the F 
core radius and exchange field. Similar oscillating behav- 
ior of T c on a ferromagnetic layer thickness has been pre- 
dicted for layered S/F structure s 12 ' 13 . Note, that the un- 
usual ground states with spontaneously formed vortices 
in multiply-connected S /F hybrids resemble the behavior 
of Josephson 7r-junctio n 14 i 15 i 16 with a step-like change in 
the F layer thicknes s 17 ' 18 ! 19 . 

It is the purpose of this paper to study the influence 



of an external magnetic field on the proximity induced 
switching between the vortex states. We focus on the be- 
havior of critical temperatures for superconducting states 
with different vorticities and, thus, we study the Little- 
Parks effect affected by the exchange interaction. 

The paper is organized as follows. In Sec. [IT] we briefly 
discuss the basic equations. In Sec. IIIII we study the 
switching between different vortex states for two model 
S/F systems placed in the external magnetic field. The 
first system consists of a thin-walled superconducting 
cylindrical shell surrounding a cylinder of a ferromag- 
netic metal. As a second example we consider a cylindri- 
cal cavity in a bulk ferromagnet covered by a thin layer 
of superconducting material. For both cases we assume 
that there is a good electrical contact between the F and 
S regions, to assure a rather strong proximity effect. We 
summarize our results in Sec. IIVI 



II. MODEL 

The calculations of the second-order superconducting 
phase transition temperature T c are based on the lin- 
earized Usadel equations 2 ^ for the averaged anomalous 
Green's functions Ff and F s for the F and S regions, re- 
spectively (see RefAi for details). The superconducting 
critical temperature T c and exchange field h is assumed 
to satisfy the dirty-limit conditions T c r <C 1 and /it < 1, 
where r is the elastic electron-scattering time. 

In the F (S) region the linearized Usadel equations take 
the form 



D f ( 2wi A 



(1) 



+ ( | w | + 1 h sgn u) F f = . 



D s ( 2iri \ , w s 

-^(^V + — Aj F s + \u,\F s =A(r). (2) 

Here Df and D s are the diffusion constants in the fer- 
romagnet and superconductor, respectively, u> = (2n + 



1) 7rT c is a Matsubara frequency at the temperature T c , 
and $o = nhc/\e \ is the flux quantum. Hereafter we 
consider only a step-like exchange field profile and, thus, 
neglect both a reduction of the magnetization in the fer- 
romagnet and magnetization leakage into the supercon- 
ductor. Such assumption is justified above T c due to the 
local nature of the exchange interaction. Calculating the 
T c value itself we can assume the superconducting order 
parameter to be vanishingly small and, thus, it is natural 
to neglect the spreading of the magnetic moment into the 
superconductor (see, e.g. ; 21 ' 22 ). 

The superconducting critical temperature T c is deter- 
mined from the self-consistency condition for the gap 
function: 



A(r) ]n^+nT c J2 

1 cO 



A(r) 



-f.(r,«) =0. (3) 



Equations fl]),© must be supplemented with the bound- 
ary condition at the outer surfaces 



d a F f>s = 0, (4) 
and at the interface between the F and S metals) 2 ^ 
a s d a F s = a f d n F f ; F 8 = F f - 7 6 £„ d a F f 



(5) 



Here £ s ( rl ) = \J ' D s ^f^/2i:T C Q is the superconducting 
(normal- metal) coherence length, at and a s are the 
normal-state conductivities of the F and S metals, 7b is 
related to the S/F boundary resistance Rb per unit area 
through 7f,£ s = RbCf, and <9 n denotes a derivative taken 
in the direction to the outer normal to the S surfaces, 
i.e., the vector n is directed from the S to the F metal at 
the S/F interface. For the sake of simplicity we assume 
h ^> ttT c0 . According to the equations |T]) , d2|) , J3]) there 
is a symmetry Ff, s (uj) = Fj s (— a->), so that we can treat 
only positive u values. 

We consider here a generic example of hybrid S /F sys- 
tems with a cylindrical symmetry: a thin-walled super- 
conducting hollow cylinder placed in an electrical contact 
with a ferromagnct (see Fig. [1]) . The ferromagnetic ma- 
terial is assumed to fill either the internal (Fig. or 
external (Fig. [TJd) region of the system. We choose cylin- 
drical coordinates (r, 9, z) as it is shown in Fig.[TJ take the 
gauge A = (0, Ag, 0), and look for a homogeneous along 
z solution of the equations ([T|) , ([2|) , (J3]) characterized by 
certain angular momentum L: 



A(r) = A(r)e 



iLO 



F f ,.(r) = f ft ,(r) 



(6) 



The vorticity parameter L just coincides with the angular 
momentum of the Cooper pair wave function. 

Certainly, the magnetization M inside the F region 
makes contribution to the vector potential A: 

B = rot A, B = H + 4ttM , 

and can modify the conditions of superconductivity nu- 
clcation in the S shell. Choosing the magnetization di- 
rection along the z— axis M = Mz , we can roughly es- 
timate the relevant change in the total magnetic flux as 
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FIG. 1: The cross section of the hybrid S/F systems un- 
der consideration: a) thin-walled superconducting shell (S) 
around a ferromagnetic cylinder (F) ; b) cylindrical cavity in a 
bulk ferromagnet (F) covered by a thin layer of superconduct- 
ing material (S). Here Rj is the radius of the F core (cavity), 
and R s is the outer (inner) radius of the S shell; (r, 8, z) is 
the cylindrical coordinate system. External magnetic field H 
applied along the z axis. 



follows: <E>m ~ An 2 Rj M. In principle, this change in the 
magnetic flux <&jy/ modifies the period of the standard 
Little-Parks oscillations and breaks the symmetry of the 
T C (H) dependence with respect to the external magnetic 
field inversion: T C (H) ^ T c (—H), However, for typi- 
cal parameters M ~ 10 2 G, T ~ 10 K, D s - 10cm 2 /s 
and Rf of order of several £/ ~ 10 nm lengths we get 
$M <C $0, and the asymmetry of the Little-Parks curve 
T C (H) may be neglected. The above effect of the magne- 
tization can be also weakened provided we decrease the 
size of the S/F system along the z— axis, going over to 
the case of a thin S/F disk when the field B is suppressed 
due to the demagnetization factor. Note that choosing 
the magnetization direction in the plane perpendicular 
to the cylinder axis we can get rid of the above magnetic 
flux correction completely because of the absence of the 
magnetization induced field component along the cylin- 
der axis. These simple estimates allow us to exclude the 
effect of magnetization on T c assuming that B ~ H and 
A e = rH/2. 

The Usadel equations (H},© take the form: 



El 

2 

Ei 
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-d r (rd r f f )-(- + -^- 
r \r 2an 

+ ihf f = 0, 

-d r {rd r fs) - f - + 

r \r Ian 

+ wf s = A, 



(7) 



(8) 



where = Q /2irH is the magnetic length. An ap- 
propriate self-consistency equation ([3]) can be rewritten 
as follows: 

A In Iz- + 2nT c V ( - - Re /.(w) ) = 0. (9) 
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III. VORTEX STATES IN THIN- WALLED 
SUPERCONDUCTING CYLINDER 

Now we proceed with calculations of the critical tem- 
perature dependence on the external magnetic field for 
two examples of multiply connected mesoscopic hybrid 
S/F systems. The first one is a ferromagnetic cylindri- 
cal filament (core) surrounded by a thin-walled supercon- 
ducting shell (see Fig.QJ,). The second one is a cylindrical 
cavity in a bulk ferromagnet covered by a thin layer of 
superconducting material (see Fig. [TJd). External mag- 
netic field H is assumed to be parallel to the cylinder 
axis (H = Hzq) for the both cases. 
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A. Ferromagnetic filament covered by a cylindrical 
superconducting shell 

Consider a superconducting cylindrical shell of a thick- 
ness W = R s — Rf -C Rf surrounding a cylinder (core) 
of a ferromagnetic metal. Here Rf is the radius of the 
F core, and R s is the outer radius of the S shell (see 
Fig. [TJi) . Naturally, to observe a pronounced influence of 
the proximity effect on the transition temperature, the 
thickness of the S shell W must be smaller than the su- 
perconducting coherence length £ s . 

The solution of Eqn. ([SJ in the F cylinder can be ex- 
pressed via the confluent hypergeometric function of the 
first kind (Kummer's function) F(a, b, z)^ 

f f (r)=Ce~^ 2 ^ L ^ 2 F(a Ll b L , 0) , (10) 

where <p is the flux of the external magnetic field H 
threading the circle of radius r in the units of the flux 
quantum $ 

cf> = 2irrAg/1>o = r 2 /2a 2 H , 

and 

2 



(11) 



a L = 



L +L + 1 



oh 




+ i ^ ; b L = \L\ + l 



(12) 



Here £/ = y/Df/h is the characteristic length scale of the 
order parameter variation in the F metal. In the dirty 
limit, the parameter £j determines both the length scale 
of oscillations and the decay length for the Cooper pair 
wave function in a ferromagnet. The boundary condi- 
tions (g]) and ^ for Eq. © take the form 



dfs_ 
dr 



0. 



(13) 



ttR 2 H/Q is the flux of the external mag- 



where x 

netic field enclosed in the F cylinder in the units of flux 
quantum $o> and 



kl(0/) = \L\ -4>f 

a L F(a L + l, b L + l, 4> f ) 



+2<t> f 



b L F(a L , b L ,(j) f ) 



(14a) 
(14b) 
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FIG. 2: (Color online) The typical dependences of the criti- 
cal temperature T c on the F core radius Rf in the absence of 
the external magnetic field for different values of the interface 
resistance 7;,. The numbers near the curves denote the cor- 
responding values of vorticity L. We choose the parameters 
(a) W = 0.5C S , &/£/ = 0.1, WO = 4.0, a s /a f = 1; (b) 
W = 0.5£ s , = 0.02, £«/£/ = 4.0, a s /a f = 0.1. 



For W <C £ s , the variations of the functions f s (r) and 
A(r) in the superconducting shell are small: / s (r) ~ /, 
A(r) ~ A. Therefore, we can average Eq. © over the 
thickness of the S shell, using the boundary conditions 
(rj"3"|) to integrate the term d r (rd r f s ). Finally, we obtain 
the following expression: 



f 



2 



R 



f 



w 



(15) 



Substituting the solution (| 1 5|) into Eq.© one obtains a 
self-consistency equation for the critical temperature Tl 
of the state with a vorticity L: 



In 



i.-O 



(16) 



where 'J is the digamma function. The depairing param- 
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FIG. 3: (Color online) The typical dependences of the critical temperature T c on the external magnetic field H for different 
values of the interface resistance 75: 75 = (o); 7;, = 0.2 (A). The magnetic field H is measured in the units of the magnetic 
flux </>/ enclosed in F cylinder. The numbers near the curves denote the corresponding values of vorticity L. Here we choose 
W = 0.5£ s ; = 0.1; = 4.0; a s /a f = 1, and different values of the F cylinder radius = (a) 0.5, (b) 1, (c) 2, 

(d) 4. The inset in panel (d) gives the zoomed part of the T C (H) line, marked by the shaded box. The dashed lines in panels 
(c, d) are guides for eye which connect the points corresponding to the T c values found for cj>f = —L, when the orbital effect in 
the depairing parameter (|17p is cancelled. 



ctcr of the mode L 



parameter kl in a simplified form: 



L + 4>f 
R 



f 



w 



(17) 



is responsible for the superconductivity destruction in the 
shell in the applied magnetic field H, due to both the 
orbital and exchange effects. As usual, the critical tem- 
perature T c of a superconductivity nuclcation in the shell 
is determined by the maximal value T^: 



T c = max{T L } 



(18) 



We start our numerical analysis from the case of zero ex- 
ternal magnetic field focusing on the effect of the S/F in- 
terface resistance Rb on the behavior of T c (Rf). Taking 
into account the asymptotic expressions for the Rum- 
mer's function F(a 7 b, z/a) for \a\ — ► 00 we obtain the 



K L (0) = \L\+u f 



h 



\L\ + 1 



(«/) 



I\L\(u f ) 



Uf 



^(1 



(19) 



The states with angular momenta ±L are degenerated 
for H = and have the same critical temperature T c . 

In Fig. [2] we present examples of dependencies of the 
critical temperature T c on the F cylinder radius Rf for 
different values of the S/F interface transparency. We see 
that for a small F cylinder radius Rf only the state 

with L = appears to be energetically favorable. The 
influence of the proximity effect is weak and the critical 
temperature T c is close to T c q. For a vortex state with 
L > 1 the T c value is suppressed because of a large or- 
bital effect. The increase in the radius Rf results in a 
decrease in T c for the state with L = and reduce the 
kinetic energy of supercurrents for L > 1. At the same 
time, the dampcd-oscillatory behavior of the supercon- 
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ducting order parameter in a fcrromagnct becomes im- 
portant, since the diameter of the F cylinder becomes 
comparable with the period of the order parameter oscil- 
lations (~ As a result, for Rf > the vortex state 
with L = 1 becomes more energetically favorable due to 
the exchange interaction. It is interesting to note that at 
small F cylinder radius (Rf < £/) the critical tempera- 
ture of the vortex free mode with L = decreases with 
an increase of the interface barrier 7b, if the condition 
IbinRf < £f is fulfilled. Such counterintiutive behavior 
is explained by the enchancement of the pair-breaking 
role of the exchange field due to the increase in the time 
of the Cooper pair stay in the F metal in the presence of 
the interface barrier—. 

Figure [3] shows examples of dependences of the crit- 
ical temperature T c on the external magnetic field H, 
obtained from Eqs. (fl~6)) . (fl~7|) , (fl"8]) for different values 
of F cylinder radius Rf. Due to the symmetry of the 
phase boundary T c (—H) = T C (H), we present here the 
curves for positive values of the external field H only. 
The phase boundary exhibits Little-Parks oscillations, 
indicating transitions between the states with different 
angular momenta L — > L ± 1 of the superconducting or- 
der parameter. For a small F cylinder radius (Fig. [3^,b) 
the influence of the exchange interaction is weak and 
the T C (H) phase boundary exhibits undamped quasiperi- 
odic oscillations as function of magnetic field and re- 
sembles the T C (H) curve for a mesoscopic disk-shaped 
superconductor^. It means that the S ring induces su- 
perconductivity in the thin F filament due to proximity, 
and the behavior of the S /F hybrid system under consid- 
eration is similar to the one for a superconducting cylin- 
der (or disk) with an inhomogeneous order parameter. 

With an increase in the F cylinder radius Rf one can 
clearly observe a shift of the main T c maximum towards 
nonzero H values. To explain this shift we note that the 
highest critical temperature for a given magnetic flux cj>f 
corresponds to the states with the angular momentum L 
close to the integer part of the —<pf value. Exactly at the 
points 4>f — —L the orbital term in the depairing param- 
eter (fTT)) is cancelled and the dependendence of T c vs tj>f 
at these points is determined only by the exchange effect. 
On the other hand, for Rf > this exchange part of the 
depairing parameter is minimal and the corresponding 
T c is maximal for a state with a nonzero vorticity (see 
Fig- HI)- Thus, the main T c maximum shifts to a certain 
nonzero <pf value. To support this explanation we have 
performed simulations of the critical temperature T c for 
the points 4>f = —L, i.e. when the orbital effect is can- 
celled. The results are shown in panels (c, d) of Fig.[3]by 
dashed lines. These curves describe a large scale (as com- 
pared to <pf) behavior of the phase boundary T C (H) and 
clearly demonstrate the shift of the main T c maximum. 
With the increase in the F cylinder radius Rf and the 
interface barrier 71, this main T c maximum shift appears 
to increase. 

Note, that a similar shift of the main T c maximum 
has been observed in mesoscopic S/F hybrid structures 



without the proximity effect^ and in the quasi-two- 
dimensional organic conductor s 2 ^ 26 . In the latter case, 
an applied magnetic field compensates the exchange 
fields of the paramagnetic ions and neutralizes the de- 
structive action of these fields 2 ^. 



B. Cylindrical cavity in a bulk ferromagnet 
covered by a superconductor 

As a second example we consider a superconducting 
shell of a thickness W = Rf — R s covering a cylindrical 
cavity of the radius Rf in a ferromagnetic material (see 
Fig. [TJa). For simplicity, we consider only rather thin S 
shells with W < £ s which allows us to assume the varia- 
tions of the functions f s (r) and A(r) in the superconduc- 
tor to be small. The solution of Eqn. ([5]) in F metal can 
be expressed via the confluent hypergeometric function 
of the second kind U(a, b, z)^, 



f f (r)=Cc-^ 2 ^ L ^ 2 U(a L , b L , 0) 



(20) 



where the flux <j>, parameters ql and bi are determined by 
the expressions (fill and (fT2l . The boundary conditions 
flU and © for Eq. © take the form 



dr 



0, * 

dr 



■X<t>f) f.(Rf), (21) 



where <f>f = irRj-H/Qo is the flux of the external mag- 
netic field enclosed in the cavity in the units of the flux 
quantum <j> , and 



Ah) 



7&£n - Rf/K L ((j)f) 

k L (4> f ) = \L\ -cj>f 

a L U(a L + l, b L + l, 4> f ) 



U(cll, b Ll <t>f) 



(22a) 
(22b) 



Assuming that the variations of the functions f s (r) and 
A(r) in the superconducting shell are small (f s (r) — /, 
A(r) ~ A) we can average Eq. © over the thickness of 
the S shell, using the boundary conditions (|2"Tj) to inte- 
grate the term d r (rd r f s ). Substitution of the solution 



f 



2 



R f 



Ah) 



w 



(23) 



into the self-consistency equation ([9]) results in the equa- 
tion for the critical temperature Tl of the state with a 
vorticity L (|15|) . where the depairing parameter Vt^ of 
the mode L is determined by the following expression: 



IT, 



cO C 2 



2 T, 



Ri 



X4>s) 



w 



(24) 



Figure 0] shows examples of dependencies of the crit- 
ical temperature T c on the external magnetic field H, 
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FIG. 4: (Color online) The typical dependencies of the critical temperature T c on the external magnetic field H for two values 
of the cylindrical cavity radius Rj ( Rf/£f = 2 (panels (a, c); -R//£/ = 4 (panels b, d) ) and different values of the the interface 
resistance 75. The magnetic field H is measured in the units of the magnetic flux <j)f enclosed in the cavity. Here we choose 
W — 0.5£ s ; = 4.0. The panels (a, b) show the case of the S/F system with a rather small ratio = 0.1 and cr s /af = 2. 

The panels (c, d) show the case of the S/F system with f s /£/ = 1 and cr 3 /cr/ = 40 . The dashed lines are guides for eye which 
connect the points corresponding to the T c values found for cj>f = —L, when the orbital effect in the depairing parameter (|24[) 
is cancelled. 



obtained from Eqs. (T5|l. 1(22]). pij) for two values of 
the cylindrical cavity radius Rf. Similarly to the previ- 
ous subsection this phase boundary resembles the T C (H) 
curve for the Little-Parks oscillations. Slow envelopes of 
the oscillating curve clearly demonstrate the shift of the 
main T c maximum towards nonzero H values. 



IV. SUMMARY 

To sum up, we have analysed the behavior of the Little- 
Parks oscillations of the critical temperature T c on an 
external magnetic field H in multiply connected S /F sys- 
tems affected by the exchange interaction. As an exam- 
ple, we have considered mesoscopic thin- walled supercon- 
ducting cylindrical shell placed in electrical contact with 
a ferromagnet. The phase-transition line T C (H) and or- 
der parameter structure have been studied on the basis 
of linearized Usadel equations. We have demonstrated 



that the exchange field provokes the switching between 
the superconducting states with different vorticities: the 
interplay between the oscillations of T c due to the Little- 
Parks effect and the oscillations due to the exchange field 
results in breaking of the strict periodicity of the T C (H) 
dependence. We have also observed a slow modulation 
of the amplitude of the quasipcriodic T C (H) oscillations. 
With an increase in the superconducting shell radius the 
envelope of the oscillating phase transition line T C {H) ex- 
hibits a shift of the main T c maximum to finite external 
magnetic field values. This shift is explained by the in- 
crease in the critical temperature of the superconducting 
states with nonzero vorticities due to the exchange in- 
teraction. The above effects strongly depend on the S/F 
interface barrier strength. We have shown that a decrease 
in the S/F transparency may stimulate the states with 
higher vorticities. 
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